Short Time Behavior in De Gennes' Reptation Model 
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To establish a standard for the distinction of reptation from other modes of polymer diffusion, we 
analytically and numerically study the displacement of the central bead of a chain diffusing through 
an ordered obstacle array for times t < 0{N^). Our theory and simulations agree quantitatively and 
show that the second moment approaches the i^' ** often viewed as signature of reptation only after a 
very long transient and only for long chains (A*' > 100). Our analytically solvable model furthermore 
predicts a very short transient for the fourth moment. This is verified by computer experiment. 
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Since the reptation model was formulated 25 years ago 
by De Gennes |1|] , it has become a widely used concept of 
polymer motion. In its original form it is concerned with 
the diffusion of a mobile polymer chain through a net of 
impenetrable and immobile obstacles. It describes this 
motion as being confined to a tube created by the sur- 
rounding obstacles. Diffusion takes place by the motion 
of little wiggles of surplus length, the "defects" , along the 
tube. Asymptotically this scenario yields simple power 
laws for various quantities M. For instance, the diffusion 
constant D, which governs the motion of the total chain 
at large time scales, should behave as D ex N~'^, where N 
is the chain length, i.e., the polymerization index of the 
macromolecule. Another important prediction concerns 
the motion of the central bead r7v/2 (t) ■ For intermediate 
times it is predicted to diffuse a mean square distance 



9i{t)^{{rN/2{t)-rN/2{0)f) 

ii/4 for To < t < T2 
ii/2 for T2 < t < T3 



(1) 



Here To = 0{N°) is a microscopic time and T2 = ©(TV^) 
is the Rouse time, i.e., the characteristic scale over which 
an unconstrained chain of length N equilibrates. The 
"reptation time" T3 = 0{N^) marks the onset of the 
asymptotic diffusional regime, where the central bead 
just follows the diffusion of the center of mass: gi{t) = 
Dt. For comparison, the Rouse model for a free nonin- 
teracting chain yields gi{t) ex t^" for t <C T2 and asymp- 
totic diffusion of both central bead and center of mass 
with diffusion constant D oc N^^ for 1^X2- 

From the beginning, the reptation model has also been 
applied to polymer motion through melts, dense solutions 
or gels, even though in such systems the surrounding 
medium more or less can relax. A lot of experimental 
or computer experimental work has aimed at verifying 
the reptation predictions in such systems. (See B for 
a recent review.) The outcome of these efforts today is 
somewhat ambiguous. The expected power laws never 



have been established beyond doubt. Compared to the 
free motion of a chain one typically finds some slowing 
down giving rise to effective power laws, with no obvious 
unique and simple interpretation pJj]. gi{t), for instance, 
typically behaves as 

5i(i)(xt^ with 0.25 < a; < 0.4 for To < i < T2 . 

It never is clear whether such deviations from predic- 
tions like (1) are due to the relaxation of the surrounding 
medium or are intrinsic to the reptation scenario in the 
available range of chain lengths. Furthermore it must 
be stressed that contrary to statements often found in 
the literature, also in a frozen disordered environment 
the asymptotic reptation behavior has not been prop- 
erly observed. In fact, recent work |g|-0 suggests that 
even strictly immobile obstacles can ruin the reptational 
power laws, due to their disordered spatial distribution. 
This finding is in contrast to the view prevailing in the 
older literature. (See g, for instance.) 

In view of this situation it is somewhat surprising to 
note that little effort has been made to clarify the im- 
plications of the original reptation model, considering 
polymer diffusion in an ordered array of fixed obstacles. 
Clearly only such work can provide the basis for a con- 
trolled analysis of effects of disorder or relaxation of the 
surrounding medium. We are aware of only two such 
studies. In early work ||^ Evans and Edwards claimed 
to find the behaviour (1), but the statistics of their data 
is insufficient for a convincing analysis. In fact, our ex- 
tended study of their model as presented below clearly 
demonstrates that for their chain lengths the t-^^^-law is 
not yet attained. More recently Deutsch and Madden 
Q reconsidered the model and found D ex N~^-^ for 
N ^ 100, in contrast to the expected iV^^-law. 

In order to establish the predictions of the pure rep- 
tation model, that serves as a starting point for the in- 
vestigation of more complicated environments, we here 
present results of an extended study of reptation through 



an ordered array of fixed obstacles in three dimensions. 
Concentrating on the internal motion in a time regime 
prior to free diffusion we present new extensive simula- 
tions, and we introduce an analytically solvable model 
that modifies the original reptation model of De Gennes 
PI to take the discreteness of the chain into account. We 
find quantitative agreement among our theory and our 
data. Simple power laws, however, are found only af- 
ter a surprisingly long initial transient and are reached 
only by quite long chains. Our theory explains the long 
transient quantitatively as a consequence of the finite seg- 
ment size. Furthermore our theory suggests to determine 
another quantity, not measured previously, which should 
show the reptational power laws more clearly, without 
the long transient. Also this prediction is verified by 
computer experiment. In summary, we provide an ana- 
lytically understood and numerically tested standard for 
reptation of discrete chains of finite length, relevant for 
the interpretation of all previous simulation results, and 
merging asymptotically {N — > oo, i ^ 1) with the ear- 
lier predicted power laws (1) ||l|. In future work, this 
should allow for a meaningful analysis of the influence 
of disorder or of some slow relaxation of the medium on 
reptation. 
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FIG. 1. The Evans-Edwards model: a) Section of a ran- 
dom walk on a square lattice, showing one hairpin. The 
crosses represent the obstacle lattice. The oval identifies 
a bead of the chain. While the hairpin diffuses past it, 
the bead is transported two steps along the tube towards 
configuration d). b) and c) give intermediate configura- 
tions. For clearer representation we have opened the base 
of the hairpins. 

In our simulations we use the Evans-Edwards model 
M: The chain is confined to a cubic lattice. A second 
identical lattice is taken to represent the obstacles and is 
placed such that the lattice points coincide with the cen- 
ters of the cells of the first lattice. Then the only moves 
allowed are those of "hairpins" , which can rotate among 
lattice bonds (see Fig. 1). This model incorporates the 
smallest possible tube diameter. It therefore shows the 
strongest restriction of chain motion and should exhibit 
reptation in clearest form. In our numerical experiment 
we concentrated on the short-time behaviour (I), search- 



ing for the t-'^/^-law. We note, however, that for shorter 
chains {N < 100) we also covered times t > T^ and re- 
covered the results of 0. A full account of our results 
will be published elsewhere. 

Fig. 2 shows gi{t) for times t ^T2. Obviously, the ex- 
pected power law gi{t) oc i^/^ indeed is attained, but only 
very slowly. It fully is developed only for surprisingly 
long chains {N > 100). It seems that for N -^ oo, there 
is a well defined limiting curve, that merges with the ex- 
pected power law only for Monte Carlo time tMC ^ 3T0*. 
(tMC = 1 stands for one attempted move per bead of the 
chain.) For our longest chain (N = 600) the i^^^'^-power 
law regime then extends over about two decades. For 
N < 100 the influence of the hairpins diffusing in from 
the ends of the chain becomes important, i.e., T2 is ex- 
ceeded, before the power law properly is developed. The 
curves then bend upwards according to Eq. (|l|). Effec- 
tive exponents j < x ^ 0.4 could easily be extracted 
for shorter chains, in full agreement with earlier observa- 
tions. 
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FIG. 2. Doubly logarithmic plot of Monte Carlo data for 
various chain lengths A''. The straight lines give slope 
1/4 or 1/3, respectively. 

To properly understand these results we worked out 
the quantitative predictions of the original reptation 
model M in a fully discretized version. The defects 
are modelled as noninteracting particles, freely diffusing 
along a lattice chain. At the chain ends defects can be 
destroyed or created, which is modelled by coupling the 
ends to large reservoirs. If a defect passes a bead, that 
bead is displaced by two steps along the "tube" , that is 
defined by the i = configuration of the chain (see Fig. 
1). Let n±(t,j) be the number of defects passing bead 
j within time interval t from the left (+) or right (— ), 
respectively. Then 2 • n{t,j) with 



nit,]) = n+{t,j) -n^{t,j) 



(2) 



gives the displacement of that bead along the "tube" . If t 
is short compared to the tube renewal or reptation time 
T3 = 0{N^), this effective diffusion for a central bead 



takes place in the fixed initial tube. Since this tube is 
a realization of a random walk on a cubic lattice with 
lattice constant a = 1, i.e., with lattice vectors s^ with 
cartesian components sf € {0, ±1} and (s^ • Sii)t = Sai, 
the average over all tube configurations yields for t -^T^ 
(a lower index c denotes the average over tube configu- 
rations and d over thermal defect dynamics) 



<?l(t) = /(rAr/2(t)-rAr/2(0))' 
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(3) 



a result that should hold up to a microstructure depen- 
dent additive constant cq much less than one. Such cor- 
rections for instance arise from configurations where the 
bead at times or i sits on the tip of a hairpin, such 
configurations being not included faithfully in our defect 
model. Note that we use the notion of the "tube" in a 
somewhat loose way, not distinguishing it properly from 
the chain configuration. In the MC-model the tube can 
be defined as the chain configuration with all hairpins 
cut off, resulting in a non-reversal random walk. The 
simple random-walk-like spatial embedding r of the in- 
ternal chain coordinates implied by Eq. (0) is justified by 
the observation that both defects and bead move along 
the chain, not along the tube. 

Our model is similar in spirit to the repton model pro- 
posed by Rubinstein [Q, but somewhat closer to the 
original De Genncs model. It allows to analytically cal- 
culate dynamical properties like the time dependent dis- 
tribution function 7'(<;n,j) of n{t,j), while taking the 
proper defect dynamics at the chain ends into account. 
Leaving details for a forthcoming paper, we only present 
some results for the central bead j = N/2. The second 
moment of the distribution of n(t, N/2) is found as 



n^{t, 
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(4) 
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where po is the density of defects. In Eq. (m we recognize 
the sum over the Rouse modes of the chain, typical for 
such problems. gi{t) [Eq. (g)] is found as a function of 
(n^) as 



9i{t) = (|2n|) 



(2(n2)) 



1/2 



2V^ + r(--,2(n2>) (5) 



2{n') 



dxx-^/^e-'' 
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where T{x,y) denotes the incomplete Gamma function. 
In Eq. (j|) both the term r(. . .) and the integral are due 



to the discrete structure of the chain. Evaluating these 
results in the time regime Tq <C i ^ T2 , that mathemat- 
ically is defined by the limit N ^ 00 and 1 ^ t < 00 
fixed, we find the intermediate asymptotics 



5l,ooW=4^"3/4^1/2^1/4 



(6) 



In the time regime r2 <C i <C T3, defined as limit iV — > 00, 
1 <C t/N^ < cx) fixed, we recover the i^'^-behavior: 
c,^(i) = 4 (po/tt)^/^ (i/7V)i/2. These are the well known 
asymptotic reptation predictions (1). 

In comparing our quantitative results to the simula- 
tion data, we have to specify the density of defects po- 
Due to the occurence of larger side loops (double hair- 
pins etc.), the mapping of hairpins onto defects is not 
one-to-one, so that pg is some effective model parameter. 
From the equilibrium statistics of a random walk it is 
found to be fixed within very close bounds: h ^ Po ^ j- 
The shape of our theoretical curves is not sensitive to the 
precise value chosen, which mainly can be absorbed into 
the time scale. We here take po = 1/4. We also need 
to fix the additive constant cq introduced above. From 



analyzing the data for microscopic times {t 



MC < 30), 

we choose to subtract cp = 0.1 from the Monte Garlo 
data for gi (t) . The only really free parameter is the scale 
T = t/tMc^ relating the time variable t of the analyti- 
cal model to the Monte Carlo time tMc- By adjusting 
this single parameter to t = 0.01, we find the results of 
Fig. 3. Note that since we devided out the intermediate 
asymptotics (6), this plot is much more sensitive to small 
deviations than the usual doubly logarithmic representa- 
tion. In view of this we believe the agreement found for 
iV > 40 to be truly remarkable. 
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FIG. 3. gi{t)/gi,oo{t) plotted against logjQ t, t = 0.01 
tMC- Full and broken lines: theory. Points: MC data. 



We now can analyze in detail the origin of the long 
transient, that in Fig. 3 extends up to t « 10^. Accord- 
ing to Eq. (||), gi{t) depends on t only via (n^(t, 7V/2)). 
Now it is found that {n'^{t, N/2)) approaches its inter- 
mediate asymptotics {n'^(t,N/2))^/'^ = {2po)^''^{t/TT)^/^ 

for times as small as i > 3. The slow transient thus is 

1 /'? 
due to the deviations of (ln|) from {2{ii?) /tt) for small 

{n?) as exhibited in Eq. (g). These deviations reflect the 
discrete structure of the chain. A simple toy model ex- 
plains this effect: The fluctuating quantity n{t,j) takes 
only integer values and its distribution rapidly looks like 
a discretized Gaussian. The width of this Gaussian, how- 
ever, increases but slowly: {n'^{t,j))^/^ ^ t^l^. Even 
for i R:! 10'^ it shows nonzero weight essentially only for 
\n\ < 5. Thus the discrete nature of the fluctuating vari- 
able n{t,j) stays important for a very large initial time 
interval. Evaluating {\n{t,j)\) with such a discretized 
Gaussian we recover the slow transient. 

These considerations immediately suggest to eliminate 
the transient behaviour by measuring (n^) instead of 
(|n|). Now for the present model it is easily checked that 
the cubic invariant gi{t) reduces to 



9i(t) 



N/2 





(7) 
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For random-walk-type chains on other lattices or in the 
continuum the spatial embedding r of the basic dynamic 
quantity (n^) differs, but (n^) always can be expressed by 
an appropriate combination of the fourth and the second 
moment of {rpf/2{t) — r^/2(0))- Fig. 4 shows a doubly 
logarithmic plot of our data for gi as compared to gi for 
chain lengths A^ = 40 and 600. Clearly our expectation is 
born out. For gi (t) only some slight initial deviation from 
the t^''* power law can be detected, that could lead to a 
tiny underestimation of the exponent (« 0.235 in place 
of 0.25). Even for A^ = 40 we find a decent i^/^-regime, 
that for A = 600 extends over almost four decades. We 
thus conclude that our fourth moment shows asymptotic 
reptational behaviour much more clearly than the com- 
monly used second moment. 
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FIG. 4. Doubly logarithmic plot of MC-data (like in Fig. 
2) for gi{t) and gi{t) and for N = 40, 600. Full line: 
i^/^'-law. Broken line: f^/^. 



In summary, combining a simple exactly solvable 
model with computer experiments, we not only recov- 
ered the asymptotic reptation prediction of a i^' ^ regime, 
but we also quantitatively explained the shape of the 
crossover function with its long initial transient. Since 
any change of the microstructure of the model like an 
enlargement of the tube radius will only increase the ini- 
tial effects, we believe our results to be typical also for 
other models. We have demonstrated that an appropri- 
ate fourth moment measuring (n^) instead of (|n|), sup- 
presses the initial transient and thus lends itself to a much 
simpler analysis. These results are hoped to provide a 
standard for the study of deviations from reptation due 
to disorder or due to a relaxation of the tube. In all cases 
we strongly recommend to include the fourth moment of 
the displacement of the central bead into the analysis. 
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